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Feedback systems which achieve subharmonic response by use of power- 
law nonlinearities and trigonometric identities are described, and the 
problem of "modes," i.e., multiple responses to the same excitation, in some 
of these systems is discussed. It is shown, by a thorough discussion of an 
example, that suitable choices of a trigonometric identity and a loop filter 
can be made which lead to locally asymptotically stable subharmonic orbits. 
An application to FM demodulation is suggested: the subharmonic modes 
can be used to reduce the index of an input wide-index wave so that a nar- 
rower IF filter than for conventional FM suffices, as in the FM demodulator 
with feedback, but without a controlled oscillator or even a mixer. 

I. INTRODUCTION 

As is known, it is possible to design feedback systems that reduce the 
index of an FM wave by an explicit use of power-law nonlinearities based 
on certain simple trigonometric identities. Circuits which use feedback to 
achieve subharmonic operation, and which are stable under changes in 
the input frequency over at least a limited range, have been built, 
tested, and described in the literature. 1 They have the behavior pre- 
dicted by the trigonometric identities. Some of them depend on non- 
linear conversion of a signal containing harmonics of 8 into one contain- 
ing only the first harmonic, and others depend on the inverse process of 
generating harmonics. 

Our purpose is to discuss the problem of "modes" (i.e., different re- 
sponses to the same excitation) in these systems, and the problem of the 
stability of the interesting subharmonic modes. We also indicate an 
application to frequency modulation: by incorporating the principle on 
which the systems are based into a frequency modulation receiver, one 
obtains circuits with some of the properties and advantages of the FM 
demodulator with feedback 23 proposed by J. G. Chaffee. However, 
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none of the circuits suggested here contains a voltage-controlled oscil- 
lator, and some do not even contain a mixer. 

To illustrate the principles involved, let us consider the trigonometric 
identity 1 

sin 30 = 3 sin — 4 sin 3 0. 

If we set 

Tf/(x) = sgn x- | x | ,/3 

the above identity can equally well be put in the form 

^(3 sin - sin 30) = 4 1/3 sin 0. 

Thus if we view, as in Fig. 1, 4 sin as the output of the nonlinearity 
\J/(- ), and feed this back through a gain of 3-4 _I/3 to an adder whose 
input is 

- sin 36, 

we will have a feedback system which is driven by —sin 36, and which 
produces sin 6, a 3-to-l reduction of index of modulation if the angle 
36 is taken to be of the form 

36 (t) = oil -f 3?(0> <P = signal. 

It is apparent 1 that other trigonometric identities can be used in 
analogous fashion to get an n-to-1 index reduction, with n any positive 
integer ^2. For example, with 

4,(x) = sgn a- \x\ 115 

we have 

^(cos 56 + 20 cos 3 6-5 cos 6) = 16 1/5 cos 6, 

corresponding to the system of Fig. 2. 
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Fig. 1 — System using 1/3-power nonlinearity. 
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Fig. 2 — Loop based on 1/5-power nonlinearity. 

Another system, based on the identity 

4 cos 3 = cos 30 + 3 cos 0, 

is shown in Fig. 3. It generates the "compressed" signal —3 cos by 
subtracting 4 cos 3 from the wide-index signal cos 30. 

II. FILTERING 

If the principles illustrated above are to be used in a communications 
receiver, it is probably desirable to perform some filtering to remove 
undesircd components of noise or feedback signal. Thus in practice the 
feedback loop would (e.g.) include a filter which removed all components 
not in the (essential) band of sin 0(- )• For many filters, and choices of 
input phase 0(-)> presence of the filter will of course mean that the 
signal in the loop is no longer so simply related to ( • ) as it was in the 
examples above: use of the trigonometric identity to relate the loop 
signal to the input may be inexact. However, there exist filters and 
choices of ( • ) for which this does not occur. (See Ref . 5, and Section V 
herein.) In any case, if the filter passes sin 0(- ) without essential distor- 
tion the identity will remain true for practical purposes. For example, 
with 



\p (x ) = sgn x • | x 



1 1/3 



again, and the identity 

iHcos 36 + 3 cos 6) = 4 1/3 cos 

we would follow \p ( • ) with a filter that passed cos but removed out- 
of-band noise, and get a system like that of Fig. 4. 



III. THE POSSIBILITY OF SEVERAL MODES 

It has been pointed out in the literature 4 that certain frequency di- 
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Fig. 3 — System based on cubic nonlinearity in the feedback. 

viders based on regeneration and modulation are not necessarily "self- 
starting," and that sizeable starting voltages may be needed to set them 
off. Put in the language of differential equations, this means for example 
that there may be two entirely different kinds of steady-state response 
(to the same steady-state signal), one (say) oscillating with large ampli- 
tudes, and the other taking place in a region of asymptotic stability 
around a critical point, with small oscillations. The "starting voltage" 
is needed to get the system out of the region of stability. 

A. J. Giger has suggested that a similar situation will obtain in realiza- 
tions of the circuits described above, even though they contain an adder 
rather than a multiplier (modulator). Clearly, whether a realization is 
self -starting is going to depend on the circuit details as well as on the 
principles at issue here, and each case will have to be studied on its own 
merits. However, these general remarks are pertinent: 

(i) The non-self -starting frequency divider is just a special case of the 
well-known but incompletely understood phenomenon that a system 
may not have a unique asymptotic response. It is a property specific to 
the systems described so far that they depend on a fractional-power 
nonlinearity, and as a result it is possible, though not necessary, that 
they fail to have unique periodic responses to some periodic signals. 5 

(ii) The desired operation of the systems proposed above depends on 
evoking a suitable subharmonic response. It is known 5 that not all solu- 
tions need contain such components of lower frequency than the input. 

(Hi) It is also known 5 that even when a subharmonic periodic solu- 
tion exists, it is itself only unique up to certain phase shifts. Specifically, 
if we obtain a solution with a component e" at when the input contains 
only harmonics of e Tinvt } n > 2, then any translation of a solution by 
2k /n for k = 1, • • • , n is also a solution. 

For example, the steady-state response of the circuit of Fig. 4 might 
or might not contain the desired subharmonics, and thus it might not 
work as planned unless care is taken to ensure that it slips into the right 
"mode" of operation initially, and that this mode is stable under the 
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perturbations due to the presence of signal and noise. The mathematical 
analysis of such phenomena is often arduous and, in many cases, quite 
incomplete. From a practical point of view the mode problem (when it 
occurs) might best be resolved by encouraging the desired mode by 
pulsing the (tuned circuit) filter, by adding automatic gain control 
features that cut in when the amplitudes of the desired mode are low, 
or by designing the filter to have zeros of transmission at certain values 
of frequency associated with the undesired modes, such as their funda- 
mentals. 

However, while it is necessary to emphasize that the problems men- 
tioned above exist, it is also important to state that the picture is not 
all black: known methods of analysis and design suffice to ensure local 
asymptotic stability of some of the subharmonic modes described above 
An example is worked out in some detail in Section V. 

IV. BIAS 

For theoretical reasons it may be undesirable, and for practical reasons 
impossible, to use a nonlinear characteristic which has an infinite slope 
at the origin. For example the singular nature of the fractional -power 
nonlinearities at the origin may preclude a conventional closed-loop, 
open-loop stability analysis by linearization around that point. Or, if 
the nonlinearities are being obtained by the use of diodes, such a slope 
is physically unattainable. These difficulties, along with most passages 
of the system through or near such a singular point, can be avoided by 
the addition of what amounts in electrical terms to "dc bias" at various 
points of the system, in such a way as to (roughly) move the operating 
point of the system to a desired region of the nonlinear characteristic 
\f/(-). In this region \J/(-) might be of Lipschitz character, or it might be 
particularly well represented by a particular diode. Such biasing can 
also be used to eliminate some equilibrium points of the system, and 
thus to reduce the number of solutions; it can also be used to increase 



cos 30 



O 



;4 3 COS 



FILTER 



l-'ig. 1 Filtering in a loop using 1/3-power nonlinearil y. 
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the size of certain regions of asymptotic stability, and thereby enhance 
the dynamic stability of desired solutions. It therefore furnishes sub- 
stantial latitude for design. 
Thus, e.g., to use the identity 

. , . 1 - cos 20 
sm 2 e = 2 ' 

we can multiply by b 2 y^ and add 

a -f- 2ab sin 

to both sides, so that 

(0 + 6 sin df = a + (6 2 /2) + 2ab sin - (b 2 /2) cos 26. 

Choosing \a\ > \b\ ensures that the right-hand side is bounded away 
from zero. Taking the square root of both sides, we base design on the 
"biased" identity 

a + b sin = (a + \b 2 + 2ab sin - \tf cos 2d) 1 ' 2 . (1 ) 

Discussion of an example of design based on (1) follows. 

V. EXAMPLE, WITH ANALYSIS OF ORBITAL ASYMPTOTIC STABILITY 

We now consider the circuit depicted in Fig. 5, with 
4>{x) = sgn x- \x\ , 
with input 

y (t) = a 2 + W - \b 2 cos 2/, 

and a filter whose impulse response k{- ) is integrable and has a Fourier 
transform 

KM = (2tt)- 1/2 f k(t)e- M dt 

J — 00 

such that K(0) = and K(l) = (27r) -1/2 (2a). For our example we 
use the second-order filter 

km- , ( yr 2a t ■ ») 

C(Za>) 2 + l(ti + C 

with c > 0, so that poles of K(- ) are in the left half-plane. 

It is now easy to show, using the theory of Fourier series, that the 
system of Fig. 5 has a subharmonic response 
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Fig. 5 — Example fur orbital asymptotic stability. 

x (t) = a + |b 2 + 2afc sin t - \b 2 cos 2«, (3) 

where x(-) is the input to the nonlincarity. The loop equation for this 
solution (with no transients!) is simply 



x(t) = y(t) + / Ht — u) sgn x(u) | x(u, 



|l/2 



dw. 



Substituting the expression that (3) gives for x(-), using the iden- 
tity (1), shows that the filter removes the dc term a from the output 
of the nonlinearity and multiplies the amplitudes of frequencies ±1 
by 2a. We remark that the constant c > in (2) does not occur in the 
subharmonic solution (3), and that the constant b does not occur in the 
transfer function K(u>). Also, for the particular input we chose, the 
presence of the filter does not render the trigonometric identity being 
used inexact. 

As is well known, for fixed values of a and c there are many ways of 
describing the circuit of Fig. 5 by differential equations so as to give rise 
to (2) as the transfer function of the filter. For some of these ways the 
solution x ( • ) found above may he orbitally asymptotically stable, for 
others it may not. In short, the stability of this subharmonic solution 
probably depends on the way in which the action of the transfer func- 
tion is represented by differential equations. We shall consider repre- 
sentations of the form 



i = Az + Sftfz + yit)), 



(4) 



where z ( • ) is a 2-vector valued function, is a 2- vector, 5 is a 2- vector 
proportional to the unit 2-vector, and A is a stable 2X2 matrix. In 
this representation the periodic solution x(-) of (3) will have the form 

x(0 = lS'z(t) +»(*). 

These representations and their simple properties are used merely as an 
illustration, because they readily admit an analysis of the asymptotic 
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stability of the orbit corresponding to the solution x ( • ) . They do not 
come near exhausting the possibilities for finding stable subharmonic 
forced oscillations in feedback systems. 
Since for t% ^ h , (4) gives 

P'z(ti) = 0'[exp A(h - h)]z(h) 

+ / 2 /3'[exp A(t 2 - u)]ty{p'z(u) + y(u))du 
J ti 

it is clear that (4 ) will represent the circuit if and only if 

/3W - A)- 1 * = (2t) ii *K(u>) = *?* 

c(ioj) 2 + *oj + c 

i.e., if and only if both 

det (iwl - A) = (t'co) 2 + c _1 ?'co + 1, 



5i — a n 



+ & 



la) — an 



— 021 



= 2aico/c. 



(5) 



^2 *W — fl 22 

The first condition is met if 

«n«22 — «i2fl2i ( = det 4 ) = 1 , 

an + «22 = — c~\ 
The second condition is equivalent to /3'5 = 2ac _1 and /3'^4 — x 5 



= 0, 



taken together. Since, to facilitate stability analysis for this example, 
we wish to impose 5i = 5 2 5* 0, the latter condition is §lA~ l l = 0, 1 = 
unit 2-vector. All these conditions together can be met in many ways. 
A convenient choice is /3 = (2a/c,0) , 8 = 1, a 2i = — a 2 i , (an — 012)022 = 
1, and 

«n(— «n — c _1 ) — a«(— On — c ) = 1, 

the last imposing a rational relation between a n and a i2 , thus leaving 
one parameter still free. 

It can be verified that when the conditions (5) above are met, then 
$ z (t ) = 2a6 sin t does define a periodic subharmonic orbit of the differ- 
ential system (4). We are now in a position to make a linear local 
asymptotic stability analysis for the subharmonic orbit (3 z(t) = 2ab sin t, 
by Lyapunov's classical theorem. 6 The (periodic) linearization matrix is 

A + diag 5/(0 



/here 
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f(l) = yf/ ( a - + 2ab sin t + \tf - \tf cos 2t) 

1 
2f (*(*)) 

If diag S = dl, d & scalar, it is easily verified that the fundamental 
matrix associated with this periodic matrix is 

*(/) = e At u(t) 

where u ( • ) > is defined by 

Uo = 1 



Un+iU) = d / f(v)u„(v)dx 

oo 

«co = S«»(o. 



We are therefore interested in the characteristic values of $>(2ir), i.e., 

e lTA u(2r). 
These are of the form 

M (2ir)e 2x " 

where /z is a characteristic value of A, and so we may conclude at once 

that if 

-Re n = — > — log w(2tt), 
2c 2-ir 

then the orbit determined by fi'z(t) = 2ab sin t is locally asymptotically 
stable, i.e., there exists a neighborhood of it from which all solutions 
approach the orbit. 
Clearly 

u(t) ^ exp {t sup | d | ■ |/(w) | ! 

0<u<t 



whence 



±- log u(2t) ^ I ' '' 



2tt & _ 2 a - | 6 

so that the inequality 
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a > c I d | + | b | 

suffices for local asymptotic stability of the orbit. 

Note that changing the sign of 6 corresponds to changing the sign in 
all the odd components in the periodic solution x(t) found above; the 
resulting function is also a periodic solution differing from x(t) only in 
phase, by exactly tt. 

VI. APPLICATION TO AN FM DEMODULATOR WITH FEEDBACK 

The examples of index reduction described in the preceding sections 
suggest that it is possible to design FM demodulators with frequency 
feedback that contain no voltage-controlled oscillators and even no 
mixers. Several methods of realizing such a possibility will be discussed, 
based on the principles exemplified. 

The simplest demodulation scheme of this sort is obtained by a 
specialization of the system of Fig. 4. We merely specify that the filter 
have a narrow passband centered around an intermediate frequency 
o) that is | the carrier frequency, and that it introduce negligible ampli- 
tude variations for signals in the passband, so as not to interfere unduly 
with the trigonometric identity. This gives rise to the demodulator of 
Fig. 6, in which 6 ( • ) has the form 

e(t) = cot + <p(t), 

with <p a baseband signal. The feedback, made at IF, reduces the modula- 
tion index 3-to-l. Here the carrier frequency is three times the inter- 
mediate frequency, but this relationship can easily be changed by re- 
modulating or, for that matter, by using a different trigonometric identity 
as the basis of design. We note that no mixer, and no voltage-controlled 
oscillator, is used. Also, the phase of the signal fed back is crucial: 
excessive phase shift in the filter is as intolerable here as in conventional 
FM with feedback. 

As a demodulator, the circuit of Fig. 6 shares with Chaffee's circuit 
the advantage that the wide band of noise which must be passed by the 
initial amplifier along with the wide-index signal is not admitted to the 
detector. This circumstance is important, because a principal object of 
feedback (in FM with feedback) is to reduce the noise level at the de- 
tector by filtering all but a small part of the noise. However, it remains 
to be seen how well the nonlinearity sgn x- J re J 1/3 performs its function 
in the presence of the wideband noise that enters it, since the resulting 
amplitude modulation a ( • ) at the input renders the trigonometric 
identity being used here inexact. This AM due to noise might be removed 
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Fig. - Possible utilization of 1/3-power nonlinearity in FM demodulator. 

by inserting a limiter and filter to follow the IF filter, as in conventional 
systems, at the expense of incurring additional phase shift. 

A final point, due to A. J. Giger, is that, unlike Chaffee's circuit, the 
present one retains the carrier phase instead of discarding it and operat- 
ing independently of it. 

VII. CIRCUITS WITHOUT FRACTIONAL-POWER NONLINEARITIES 

It is straightforward to generate other, quite different designs based 
on the same identity, 

cos 30 = 4 cos 3 — 3 cos 0, 

designs which do not depend on a fractional-power nonlincarity, and 
so do not incur the problems above. For example, Fig. 7 shows a design 
very much like that of the conventional FM with feedback demodula- 
tor, except that the detector-controlled oscillator is replaced by the 
nonlinearity 

Aif - 3y. 

Applied to cos 0, this gives cos 30, to be used as the feedback input to 
the mixer. If the other mixer input, i.e., the incoming signal, is the 
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Fig. 7 — Cubic nonlinearity used in FM feedback system with modulator. 
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wide-index signal cos 40, the filter can be made to select the difference 
frequency component cos to complete the loop and provide the feed- 
back signal. In this system the carrier frequency is four times the inter- 
mediate frequency. 

A particularly simple circuit, based on Fig. 3 and using only an adder, 
is depicted by Fig. 8. In this design the feedback is through the simple 
cubic nonlinearity 

27 

If —3 cos is applied to this, and the output is combined in the adder 
with an incoming wide-index signal cos 30, the adder output is —3 cos 
0. This passes substantially unchanged through the narrow IF filter 
suitable for the low-index wave cos 0, while of the wide band of noise ac- 
companying cos 30 at the input only a narrow band can pass the IF filter. 
For practical purposes, the cubic characteristic would only be required 
over the range | x | ^3, and standard stability analyses can be used. 
Again, carrier phase is retained. 
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Fig. 8 Cubic nonlinearity used in FM feedback system with adder. 
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